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Abstract 

Two neutrino double p decay can create irremovable background even in high energy resolution detectors search- 
ing for neutrinoless double decay due to random coincidence of 2v2[5 events in case of poor time resolution. 
Possibilities to suppress this background in cryogenic scintillating bolometers are discussed. It is shown that the 
present bolometric detector technologies enable to control this form of background at the level required to explore 
the inverted hierarchy of the neutrino mass pattern, including the case of bolometers searching for the neutrinoless 
double p decay of "'"Mo, which is characterized by a relatively short two neutrino double p decay half-life. 



1. Introduction 

Neutrinoless double beta (0v2j8) decay is a key pro- 
cess in particle physics since it provides the only ex- 
perimentally viable possibility to test the Majorana na- 
ture of neutrino and the lepton number conservation, 
establishing in the meantime the absolute scale and the 
hierarchy of the neutrino masses [1]. 

One of the most important goal of the next gener- 
ation double p decay experiments is to explore the in- 
verted hierarchy of the neutrino masses. In the inverted 
scheme the effective Majorana {m^) neutrino mass is 
expected to be in the interval ~ 0.02 - 0.05 eV. In or- 
der to check this range, an experimental sensitivity (in 
terms of half-life) for the most promising nuclei should 
be at the level Txji ~ 10^^ - 10^^ yr. This requires a 
detector containing a large number of studied nuclei 
(~ 10^' - 10^*^), with high energy resolution (at most 
a few % at the energy of the decay Q2/i), large (ideally 
100%) detection efficiency, very low (ideally zero) ra- 
dioactive background. 

Besides HPGe detectors used to search for 0v2j6 de- 
cay of ^^Ge [2, 3], cryogenic bolometers [4, 5] - lumi- 
nescent [6, 7, 8, 9, 10] or not [11, 12, 13] - are excellent 
candidates to realize large-scale high-sensitivity exper- 
iments involving different isotopes with high energy 
resolution (a few keV) and detection efficiency (near 
70% - 90%, depending on the crystal composition and 
size). 

The CUORE cryogenic experiment [12, 13], built 
on the successful Cuoricino and searching for 0v2y6 de- 
cay of '^"Te with the help of TeOa detectors, is by far 
the most advanced bolometric search and is under con- 
struction now, while several searches (among which 
LUCIFER [14] and AMoRE [15]), aiming to use dif- 
ferent luminescent bolometers to search for double p 



decay of '^^Se (ZnSe [14, 16]), "^Cd (CdW04 [17]), 
i™Mo (CaMo04 [15] and ZnMo04 [18, 19, 20, 21]), 
and '^°Te (TeOi [22]) are in the R&D stage. 

However, a disadvantage of cryogenic bolometers is 
their poor time resolution. This can lead to a back- 
ground component at the energy Qi/? due to random 
coincidences of lower energy signals, especially those 
due to the unavoidable two neutrino double p decay 
(2v2j6) events. 

The random coincidence of 2v2y6 events as a source 
of background in high sensitivity 0v2y6 experiments 
was considered and discussed for the first time in [20]. 
In this work, the contribution of random coincidences 
(rc) of 2v2y6 events to the counting rate in the energy 
region of the expected Qvip peak is estimated. Meth- 
ods to suppress the background are discussed. 



2. Random coincidence of 2v2/3 events 

Energy spectra of p particles emitted in 2v2;6 decay 
are related to the 2-dimensional distribution pi2(fi, fa) 
(see f.e. [23] and refs. therein) 

PnitiJi) = eip\F{tuZ) ■ e2PiF{t2,Z) ■ (to - f i - fa)^ 

(1) 

where f, is the kinetic energy of the ;-th electron (all 
energies here are in units of the electron mass otqc^), Iq 
is the energy available in the 2p process, is the mo- 
mentum of the ;-th electron p, = ^fti(tp^Tj (in units 
of niQc), and e, = f, -i- 1 . The Fermi function F{t, Z), 
which takes into account the influence of the electric 
field of the nucleus on the emitted electrons, is defined 



F{t,Z) — const 



p^" ^exp{7iT]) 



T(s + irf) 



(2) 
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where s = ^Jl - {aZf-, jj = aZe/p, a = 1/137.036, Z 
is the atomic number of the daughter nucleus (Z > 
for 2/3^ and Z < for 2j6^ decay), and T the gamma 
function. 

The distribution p(t) for the sum of electron energies 
t = ti + t2 is obtained by integration 

p(f)= r Pn{t-t2,t2)dt2. (3) 
Jo 

This distribution is shown for '*'Mo in Fig. 1. 
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Figure 1 : (Color online) Distribution for the sum of energies of two 
electrons emitted in 2v2/? decay of '""Mo and energy spectrum of 
10* two randomly coincident 2v2/3 events for '""Mo obtained by 
Monte Carlo sampling. The approximation of the random coinci- 
dence spectrum by the expression (7) is shown by the solid(red) line. 



The Primakoff-Rosen (PR) approximation for the 
Fermi function F{t,Z) ~ e/p [24], which is adequate 
for Z > 0, allows to simplify Eq. (1) to the expression 

Pn(ti^t2) = (fi + lf{t2 + If (to -H- t2f (4) 

and to obtain the formula for p(f) analytically: 

p™(f) = f(fo - tfit'' + lOr^ + AQt- + 60f + 30). (5) 

The energy distribution for two randomly coincident 
2v2y6 decays prdt) can be obtained by numerical con- 
volution 

(6) 



Prcif) - pit- X)p(x) dx, 

Jo 

or with a Monte Carlo method by sampling energy re- 
leases in two independent 2v2y6 events in accordance 
with the distribution (3) and adding them. The energy 
spectrum obtained by sampling 10** coincident 2v2/3 
events for '™Mo is shown in Fig. 1. (We assume here 
an ideal energy resolution of the detector) This distri- 
bution can be approximated by the following compact 
expression,' similar to that reported in Eq. (5); 



p,,(f) = f'(2fo-f)'° 



(7) 



However, the coefficients a, are different for differ- 
ent isotopes; for '^^Se, '™Mo, "''Cd, '^°Te (which are 
the nearest aims of the bolometric 2y3 experiments). 



Table 1: Coefficients a, in the energy distribution (7) for two ran- 
domly coincident 2v2/? events for '^^Se, '"Mo, ' '*Cd. '^^"Te. 



Isotope 





**2Se 


ioomo 


""Cd 




flo 


3446.59 


5827.48 


20093.8 


15145.6 


fll 


-7746.37 


-14399.7 


-1318.65 


5554.78 


ai 


22574.7 


34128.1 


69134.6 


39930.0 


fl3 


-16189.1 


-23815.5 


-31971.3 


-13338.5 


fl4 


8467.01 


11271.7 


17976.8 


7689.45 


fl5 


-2156.83 


-2711.31 


-3486.40 


-813.887 


fl6 


337.172 


390.396 


406.762 


-80.3126 


fl7 


-28.9146 


-30.8774 


-30.5846 


19.1035 


fl8 


1. 


1. 


1. 


-1. 



they are given in Table 1 . The approximation for ""^Mo 
is shown in Fig. 1. 

The random coincidence counting rate in a cho- 
sen energy interval AE is determined by the time reso- 
lution of the detector r and the counting rate for single 
2v2y6 events Iq: 



Ire 



ln2N 



(8) 



where is the number of 2(3 decaying nuclei under in- 
vestigation, and e is the probability of registration of 
events in the AE interval. In Eq. (8) and in the follow- 
ing, we assume that if two events occur in the detector 
within a temporal distance lower than the time resolu- 
tion T they give rise to a single signal with an amplitude 
equal to the sum of the amplitudes expected for the two 
separated signals. The calculated probabilities at the 
energy Qa/j of the Ov2/3 decay for AE - 1 keV interval 
are equal to e = 3.5 x 10""^ for '^^Se and e = 3.3 x 10""^ 
forio"Mo,ii''Cd, i^^OTe. 

Counting rates of detectors with 100 cm"* volume 
(typical for large mass bolometers) at the energy of 
0v2y6 decay for different 2yS candidates and compounds 
are presented in Table 2. We assume 100% isotopical 
enrichment for ^^Se, "^°Mo, "^Cd, and natural abun- 
dance (34.08%) for '■'°Te. The time resolution of a 
detector is assumed as t = 1 ms. The reported rates 
scale linearly with the time resolution. 

Background B,c caused by random coincidences of 
2v2j0 events has the following dependence on the en- 
ergy resolution R, the volume of the detector V, the 
abundance or enrichment 6 of the candidate nuclei con- 
tained in the detector: 



Brc 



t R- (rf;f )-2 . 



(9) 



' One can expect a polynomial of 2 1 -st degree because of the for- 
mulae (5) and (6). 
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One can conclude from Table 2 that the most impor- 
tant random coincidence background is for "^"Mo due 
to the relatively short 2v2y8 half-life. However, a non- 
negligible contribution is expected also for other iso- 
topes in case of big volume and poor time resolution 
of the single detectors. 

Obviously any other source of background with en- 
ergy high enough can contribute due to random coinci- 
dences. For example, the presence of 234Hip^ (belong- 
ing to the series of ^^^U) with a relatively high activity 
of 1 mBq/kg in 100 cm"* ZnMo04 crystal will result 
in additional background due to random coincidence 
with 2v2/3 events of 3.8 x 10"^ counts/(keV-kg-y) at 
loojVIo Q213 energy. This contribution is however much 
less important than that due to the 2v2y6 decay alone. 
In addition, our simulations show that coincidence of 
2v2y6 signals with low energy events is not problem- 
atic, because of the quite steep shape of 2v2yS spec- 
trum near Qa^g. Low energy signals due to spurious 
sources, like microphonic noise, which can in princi- 
ple contribute with a high rate, are in general easily 
rejected in bolometers thanks to pulse shape discrim- 
ination. In conclusion, while a generic source can be 
reduced by careful shielding, purification of materials, 
improvement of the noise figure and anti-coincidence 
technique, the random coincidence of 2v2y6 events is 
a background hard to suppress, as it is related to the 
presence itself of the isotope under investigation. The 
possibilities to decrease this type of background at an 
acceptable level are discussed in the next Section. 

3. Pile-up rejection in scintillating bolometers 

A cryogenic bolometer [26] consists of an energy 
absorber (a single diamagnetic dielectric crystal in 
Qvlfi applications) thermally linked to a temperature 
sensor, that in some cases may be sensitive to out of 
equilibrium phonons. The heat signal, collected at very 
low temperatures (typically < 20 mK for large bolome- 
ters), consists of a temperature rise of the whole detec- 
tor determined by a nuclear event. 

The majority of the most promising high Qayj-value 
(> 2.5 MeV) candidates can be studied with the bolo- 
metric technique in the "source=detector" approach, 
joining high energy resolution and large efficiency [4, 
5]. Ultra-pure crystals up to 100 - 1000 g can be grown 
with materials containing appealing candidates. Ar- 
rays of the single crystalline modules allow achieving 
total masses of the order of 100- 1000 kg [12, 13, 14], 
necessary to explore the inverted hierarchy region. 

An excellent choice for the bolometric material, per- 
formed in the Cuoricino and CUORE experiments [11, 
12], consists of TeOi (tellurite) that has a very large 
(27% in mass) natural content of the 0v2y6 candidate 
'^°Te. In terms of background, the experience pro- 
vided by Te02 searches [27] shows clearly that energy- 
degraded a particles, emitted by the material surfaces 



facing the detectors or by the detector surfaces them- 
selves, are expected to be the dominant contribution 
in all high (22/3-value candidates, for which the sig- 
nal falls in a region practically free of y background. 
The a background component however can be made 
negligible using scintillating or in general luminescent 
(including Cherencov light [10, 22]) bolometers. In 
fact, since the a light yield is generally appreciably 
different from the Ply light yield at equal deposited en- 
ergy, while the thermal response is substantially equiv- 
alent, the simultaneous detection of light and heat sig- 
nals, and the comparison of the respective amplitudes, 
represents a powerful tool for ajfi discrimination and 
therefore for a background rejection [6, 7, 8]. Scin- 
tillation photons are usually detected by a dedicated 
bolometer, in the form of a thin slab, opaque to the 
emitted light and equipped with its own temperature 
sensor. The light absorber, normally a Ge, Si or Si- 
coated AI2O3 slab, is placed close to a flat face of the 
main scintillating crystal. 

A wealth of prelimimary experimental results [9, 
16, 17, 18, 19, 20, 21] show that this method is ef- 
fective and a rejection factors even much better than 
99.9% can be achieved. Once that this discrimination 
capability of scintillating bolometer is taken into ac- 
count, a detailed background analysis, based on rea- 
sonable assumptions on internal radioactive contam- 
ination, shows that a residual background level of ~ 
10""* counts/(keV-kg-y)canbe safely assumed [20, 21], 
opening the opportunity to explore the inverted hier- 
archy region of the neutrino mass pattern. However, 
the random coincidences of 2v2y8 events discussed in 
Section 2 need to be kept under control. According to 
Eq. (8), the rate of rc-2v2y8 events is proportional to 
the time resolution of the detector Therefore the time 
properties of the signal from a cryogenic detector play 
a crucial role in this form of background. 

The large mass (~ 800 g), high energy resolution 
(~ 3 -4 keV FWHM at 2615 keV) detectors developed 
for Cuoricino and CUORE [11, 12] represent a sort of 
paradigma for Ov2/3 decay bolometers, inspiring also 
other proposed experiments and the related R&D ac- 
tivities. As temperature sensors. Neutron Transmuta- 
tion Doped (NTD) Ge elements are used, characterized 
by a high impedance (1-100 MQ) and a high sensitivity 
{-dlogR/dlogT ~ 10). Other promising solutions for 
the thermal sensors have been used [15, 28, 29, 30], but 
for the moment only the energy resolution and the re- 
liability provided by NTD-Ge-based bolometers seem 
to be compatible with a large-scale 0v2;6 experiment. 
In these devices, the NTD Ge thermistor is glued to 
the Te02 crystal by means of a two-component epoxy. 
Their time resolution is strictly related to the thermal- 
signal risetime. 

The temporal behaviour of the thermal pulse, and 
therefore its risetime, can be understood thanks to a 
thermal model for the whole detector, described else- 
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Table 2: Counting rate of two randomly coincident 2y2y3 events in cryogenic Zn**-Se, '*''Ca""'Mo04, Zn""'Mo04, "^CdW04 and Te02 
detectors of 100 cm^^ volume. Enrichment of '^^Se, '""Mo and ' '*Cd is assumed 100%, while for Te the natural isotopic abundance (34.08%) 
is taken. C is the mass concentration of the isotope of interest, p is the density of the material (g/cm'), N is the number of 2/3 candidate nuclei 
in one detector, and B„. is the counting rate at Qip (counts/(keV-kg-yr)) under assumption of 1 ms time resolution of the detectors. 



Isotope 


rf;f (yr) [25] 


Detector (p) 


C 




Brc 






9.2x10"'' 


Zn^^Se (5.65) 


55.6% 


2.31x1024 


5.9x10- 


-6 


ioomo 


7.1x10"* 


40Ca"'0MoO4 (4.35) 


49.0% 


1.28x1024 


3.8x10- 


-4 






Zn"'°Mo04 (4.3) 


43.6% 


1.13x1024 


2.9x10- 


-4 




2.8x10''' 


i"'CdW04 (8.0) 


31.9% 


1.32x1024 


1.4x10- 


-5 


i30Te 


6.8x1020 


TeOz (5.9) 


27.2% 


0.76x1024 


1.1x10- 


-8 



where [31, 32, 33, 34]. The model predicts that large 
mass detectors with NTD-Ge readout have risetimes 
of the order of tens of milliseconds. This is confirmed 
experimentally in the Cuoricino detectors, which ex- 
hibited pulse risetimes of the order of ~ 50 ms [11]. 
Similar values are expected for any cryogenic bolome- 
ters with a volume of the order of 100 cm^ and 
based on NTD Ge thermistors. Faster risetimes could 
be observed if an important component of the en- 
ergy reaches the thermistors in the form of athermal 
phonons, but this possibility depends critically on the 
nature of the main crystal and of the crystal-glue- 
thermistor interface. The most conservative approach 
consists in assuming the slow risetime evaluated and 
observed in TeOa bolometers. 

The time resolution of the main crystal can be made 
substantially shorter than the risetime, taking advan- 
tage of the excellent signal-to-noise ratio expected at 
the 0v2y6 energy [35], which is of the order of 2000 to 1 
in Te02 crystals. (However, this high value is still to be 
proved for other bolometric materials relevant for v2y3 
decay.) Even though the time resolution proved shorter 
by a factor 10 with respect to the present TeOa risetime 
values, and therefore around 5 ms, the background 
values reported in Table 2 should be multiplied by a 
factor 5, bringing them above lO^^* counts/(keV-kg-yr) 
for 'O^Mo. Background due to random coincidence of 
2v2/3 events would be then dominant. The advantage 
of scintillating bolometers, which promise to keep the 
other sources around or below 10^4 counts/(keV-kg-yr) 
[20, 21], would be substantially compromised. 

However, the simultaneous detection of light and 
heat which characterizes scintillating bolometers offers 
the possibility to control this problematic background 
source as well. In fact, a much faster risetime in the 
light-detector signal is expected, due to the much lower 
heat capacity of the energy absorber, which has now 
a mass of only a few grams at most. Following the 
bolometer thermal model [31, 32, 33, 34], the light- 
detector risetime can be reduced down to ~ 1 ms. As- 
suming a similar time resolution, the 2v2j6 background 
contribution is in the ranges shown in Table 2 thanks 



to the fast response provided by the Ught detector. 
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Figure 2: (Color online) Simulated piled-up pulses (solid/red lines) 
using real pulse shape and noise from a working light detector cou- 
pled to a ZnMo04 scintillating bolometer: (a) pile-up of two pulses 
shifted by 3 ms with equal amplitude; (b) pile-up of two pulses 
shifted by 3 ms with amplitude ratio equal to 4 (the smaller pulse oc- 
curs first); (c) a single pulse. The typical single-signal pulse shape, 
obtained by fitting an average pulse, is plotted as well (dashed/blue 
Hues). In all cases, the signal-to-noise ratio is that expected for a 
0v2/3 signal. The difference in shape between piled-up and single 
pulses is small, especially for unequal amplitudes, but appreciable. 

The above discussion is simplified: the capability 
to discriminate two close-in-time events cannot be re- 
duced to a single parameter such as the detector time 
resolution t. In fact, it depends smoothly on the tempo- 
ral distance between the two events and on the ratio be- 
tween their two amplitudes as well. Furthermore, the 
pile-up discrimination capability is influenced by the 
signal-to-noise ratio in the light detector [35] and by 
the pulse shape and noise features. A complete analy- 
sis, whose details will be reported elsewhere in a ded- 
icated publication, has been performed. In this letter. 
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we present the main results of this investigation and the 
most important conclusions. 

An experiment based on modules of Zn"'°Mo04 
crystals was considered. We have used pulse shapes 
and noise from a real light detector, of the same type as 
those described in [19, 20], coupled to a ZnMo04 scin- 
tillating crystal. The pile-up phenomenon was studied 
by generating light pulses with the observed experi- 
mental shape on top of experimental noisy baselines. 
In particular, pair of pulses were generated with ran- 
dom time distances with a flat distribution up to 10 ms 
- in fact, the interarrival time distribution is practically 
constant over the [0, 10] ms range (this is the rele- 
vant time interval since it allows to fully explore the 
most problematic pile-up case, i.e. that occurring on 
the pulse risetime which is of the order of 3 ms). In 
Fig. 2, two pile-up emblematic cases (both with 3 ms 
time separation) extracted from the performed simula- 
tion are shown, along with a single pulse as a reference. 

As a first step, we defined a 90% efficiency in ac- 
cepting a pulse from the light detector as a potentially 
good Qv2/3 pulse using opportune signal filtering and 
three different pulse-shape indicators: (i) the risetime 
from 15% to 90% of the maximum amplitude; (ii) the 

evaluated using an average pulse as a standard shape 
function; (iii) the pulse shape parameter defined in [36] 
which also uses a standard pulse-shape function. The 
rejection efficiency of piled-up pulses was then tested. 
In each pulse pair, the amplitude of the first pulse Ai 
was extracted by sampling the 2v2y6 distribution, while 
the amplitude of the second pulse A2 was chosen as 
G2/?(""'Mo)-Ai + AE, where AE is a random compo- 
nent in the interval [-5,-1-5] keV. 

The generated pulse amplitudes were chosen so as to 
fix the signal-to-noise ratio at the level expected for a 
0v2;6 signal, i.e. of the order of 30, as shown in Fig. 2; 
in fact, the typical light energy collected by the light 
detectors in ZnMo04 scintillating bolometers realized 
so far is of the order of 1 keV for 1 MeV energy in the 
heat channel [18, 19, 20, 21], while the typical RMS 
noise of the light detector can be conservatively taken 
as 100 eV, although values as low as 30 eV were ob- 
served [19]. 

The piled-up pulses generated in the simulation 
were analyzed with the mentioned pulse-shape indica- 
tors. Using the risetime (after low-pass filtering), an 
excellent pile-up rejection efficiency was obtained. A 
comparison between the risetime distribution for gen- 
uine single pulses and piled-up pulses generated as de- 
scribed above is shown in Fig. 3. More quantitatively, 
the same procedure that retains 90% of genuine sin- 
gle pulses rejects 80%-90% of piled-up pulses when 
their sum amplitude is in the region of Q2/S and the 
diff'erence between the arrival times of the two pulses 
covers uniformly the interval [0, 10] ms. For exam- 
ple, the analysis of the sample reported in Fig. 3 ex- 
cludes 83% of piled-up pulses when accepting 90% of 



good pulses. The other two indicators provide equiva- 
lent or even better results. However, we prefer here to 
consider conservatively the results obtained with the 
method of the risetime, as this parameter is an intrin- 
sic property of each signal that does not require the 
comparison with a standard shape. This comparison 
in fact implies a delicate synchronization between the 
single pulse and the standard-shape pulse that will be 
discussed in the mentioned more complete work. 

The results of the simulation show that the contribu- 
tion to the background of the piled-up events is sub- 
stantially equivalent to that obtained when assuming a 
time resolution t of 1 ms in Eq. 9 (since 80% - 90% 
of the pulses are rejected in the region of Qv2/3 decay 
inside a pile-up relevant range of 10 ms), and therefore 
confirming the evaluation for ZnMo04 reported in Ta- 
ble 2. 

800 
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Figure 3: (Color online) Risetime distribution for two populations 
of 5000 generated events each. The soHd (red) line refers to single 
pulses; the dashed (blue) line is obtained with piled-up pulses sepa- 
rated by a time distance coveting unifomily the interval [0, 10] ms, 
with amplitudes sampUng the 2v2fi spectrum and adding so as to fall 
in the region of the 0v2/} expected peak. 

We can conclude therefore that light detectors at the 
present technological level are compatible with next- 
generation 0v2y6 decay experiments based on ZnMo04 
crystals with background in the 10""* counts/(keV-kg-y) 
scale, confirming that this class of experiments has the 
potential to explore the inverted hierarchy region of the 
neutrino mass pattern [20, 21]. 

4. Conclusions and prospects 

Random coincidence of 2v2y6 events is an irremov- 
able background source in a large scale 2/3 experiments 
using detectors with slow response time, such as large 
mass cryogenic bolometers with NTD Ge readout. 

Advancement of time resolution of cryogenic de- 
tectors plays a key role to suppress the background. 
However, we have shown that the present technology 
is already compatible with searches at the sensitivity 
frontier For further improvements, experimental ef- 
forts should be concentrated on the time properties of 
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the light signal, potentially much faster than the heat 
pulses of a scintillating bolometer. Achieving a time 
resolution below 0.1 ms could make the background 
totally negligible even for the difficult case of ""^Mo. 
Such a performance could be obtained by using sensors 
sensitive to out-of-equilibriumphonons or intrinsically 
fast [28, 29, 30]. 

A more direct way to decrease the pile-up effect is to 
reduce the volume of the main absorber (and increas- 
ing correspondingly the number of array elements), on 
which the random concindence rate depends quadrat- 
ically, as shown in Eq. (9). Cryogenic detectors with 
space resolution could allow to reduce the background 
further 
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